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[57] ABSTRACT 

A synthetic aperture radar (10) is employed for Ak 
measurement of ocean current from a spacecraft (11) 
without the need for a narrow beam and long observa- 
tion times. The SAR signal is compressed (12) to pro- 
vide image data for different sections of the chirp band- 
width, equivalent to frequencies fi(t G ,t), f 2 (L,t) . . . 
f/i(t a ,t), and a common area for the separate image fields 
is selected (14). The image for the selected area at each 
frequency is deconvolved (16) to obtain the image sig- 
nals for the different frequencies (fi, f 2 . . . f w ) and the 
same area. A product of pairs of signals is formed (18, 
20), Fourier transformed (22) and squared (24). The 
spectrum thus obtained from different areas for the 
same pair of frequencies fy, f )+ n ,k are added (26) to 
provide an improved signal to noise ratio. The shift of 
the peak from the center of the spectrum is measured 
and compared (28) to the expected shift due to the phase 
velocity of the Bragg scattering wave. Any difference is 
a measure of current velocity v c (Ak). 

4 Claims, 4 Drawing Figures 
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METHOD AND APPARATUS FOR AK SYNTHETIC 
APERTURE RADAR MEASUREMENT OF OCEAN 
CURRENT 

5 

ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 10 
(72 Stat. 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to the use of a chirped syn- 
thetic aperture radar (SAR) to measure ocean current 15 
velocity from spacecraft, and more particularly to a 
SAR system for measurement of surface currents of the 
ocean from the ocean wavelength of 27r/Ak, where Ak 
is the difference in wavenumber .between two frequen- 
cies fi and $2 of a chirped SAR. 20 

Surface currents of the ocean vary widely from the 
current velocity below about one meter due to local . 
surface wind and wave fields that differ from the nor- 
mal waves resulting from geostrophic forces and tides. 

It would be desirable to map the deeper ocean current 25 
velocity from a spacecraft not only for climate predic- 
tion in coastal areas, but also for navigation and man- 
agement of fishing operations. 

Previous work has shown that Ak radar techniques 
can measure current velocities, and determine dominant 30 
wavelengths, based on the fact that if the Bragg angle 6 
between the incoming radiation and the mean surface is 
larger than the rms surface slope, but smaller than about 
70°, microwave scattering is well defined by the Bragg 
law, as discussed by W. J. Plant and D. L. Schuler in a 35 
paper titled “Remote sensing of the sea surface using 
one- and two-frequency microwave techniques” pub- 
lished in Radio Science 15, 605 (1980). 

In that paper, the authors describe a dual-frequency 
technique. Two coherent microwave signals are trans- 40 
mitted from a stationary CW radar, and signals scat- 
tered from the same spot at the two frequencies are 
multiplied to produce a spectrum with a sharp line, the 
Ak line. The authors show that for a narrow beam and 
relatively long observation times, a dual-frequency sys- 45 
tern may thus be used to measure current velocities and 
determine dominant wavelengths. 

SUMMARY OF THE INVENTION 

An object of this invention is to use the Ak technique 50 
for determining ocean current velocities using a chirped 
SAR signal from a spacecraft radar system, without the 
need for a narrow beam and long observation times. 

In accordance with the present invention, chirped 
synthetic aperture radar signals S(t a ,t) are first com- 55 
pressed to provide signals e/t CD ,t 0 ) for different sections 
of the chirp bandwidth equivalent to signals obtained 
for different frequencies fi(t fl ,t), f 2 (L,t) . . . fi,(t*,t), 
where S(t a ,t) represent distinct chirped radar input sig- 
nals at times t a and t, i.e., where S(t^,t) is a two-dimen- 60 
sional signal. A common area, K, for the separate image . 
fields is selected to produce signals e y *(Wt 0 ), where the 
cap (*) over the e indicates the total signal from area. 
The cap, which is then carried over all symbols for 
subsequent signals carries the same meaning. Note that 65 
the upper case letter K used for area does not refer to 
the Ak Bragg wave, and that the lower case letter k in 
subscripts refers to this area K, but that Ak does refer to 
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difference in the wavenumber of an ocean wave of 
wavelength related to 27r/Ak. The common rectangular 
area, K, is of dimensions UxU centered at t ao k, t^for 
the separate image fields at the different frequencies. 
The image zjk{t ao >to) is deconvolved to obtain the signals 
S/*(Wt 0 ) for the different frequencies and the same 
area. A product of pairs of signals S jk and S is then 
produced that is a two-dimensional signal function of 
Ak and t ao , t 0 - The product is processed through a Fast 
Fourier Transform (FFT) and a squarer. The process is 
repeated for different areas of the same frequencies. The 
spectrum thus obtained for different areas using the 
same pair of frequencies are added to provide an im- 
proved signal to noise ratio for the spectrum. A Dop- 
pler shift fcof the spectrum peak from the center of the 
spectrum is compared to the expec ted Doppler shift due 
to the phase velocity v /7 =Vg/Ak of the Bragg scatter- 
ing wave, i.e., actual velocity without current where g 
is the acceleration due to gravity. This difference is due 
to the current velocity v c (Ak) as effecting the ocean 
wavelength of 27r/Ak. In other words, from the equa- 
tion 

where fa is the radar frequency and C is the speed of 
light, the total velocity V(Ak) is determined. The radial 
component of velocity v c (Ak) is then determined from 
the equation V(Ak)=v^+v c (Ak). 

The novel features that are considered characteristic 
of this invention are set forth with particularity in the 
appended claims. The invention will best be understood 
from the following description when read in connection 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a spacecraft having a chirped syn- 
thetic aperture radar (SAR) illuminating the ocean for 
the purpose of determining the ocean current from the 
Ak Bragg condition. 

FIG. 2 illustrates the diagram of FIG. 1 looking in the 
direction of the spacecraft flight path. 

FIG. 3 is a functional block diagram of a Ak SAR 
processor in accordance with the present invention. 

FIG. 4 is a mathematical flow chart for the operations 
of the processor of FIG. 1 which may be carried out in 
a digital SAR processor. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Referring to FIG. 1, the chirped side looking syn- 
thetic aperture radar 10 is carried by a spacecraft 11 
over the ocean. As shown, the input signals S(t a ) and 
S(t) appear to be vertical, but in practice will be at a 
slant angle /3 that is the complement (90° — 0) of the 
Bragg angle 0 between the incoming radiation and the 
mean surface, where 0 is larger than the rms surface 
slope, but smaller than about 70° as shown in an end 
view in FIG. 2, i.e., shown in a view looking in the 
direction of the spacecrafts flight path. The SAR signal 
is processed to provide images for different sections of 
the chirp bandwidth, equivalent to images obtained for 
different frequencies. A common area for the separate 
image fields is selected and the image deconvolved to 
obtain the signal for different frequencies and the same 
area. A product of pairs of signals is Fourier trans- 
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formed and squared, and the spectrum thus obtained for 
different images for the same pair of frequencies added 
to provide an improved signal to noise ratio. The shift 
of the peak from the center of the spectrum is deter- 
mined and compared to the expected shift due to the 5 
phase velocity of the Bragg scattering wave. Any dif- 
ference is due to the radial component of the current 
velocity as affecting the ocean wavelength, and thus a 
measure of ocean current. 

For simplicity, only the dominant waves of deeper 1° 
ocean currents due to geostrophic forces and tides are 
shown in FIGS. 1 and 2, it being understood that sur- 
face waves due to winds and other factors are also 
present and contribute to the input signal S(t fl ,t). How- 
ever, these surface waves to not interfere with measure- 15 
ment of the ocean currents by this SAR system which 
effectively obtains a narrow beam using SAR process- 
ing and a high signal to noise ratio using spatial averag- 
es- 90 

FIG. 3 illustrates the SAR processing system of the 
present invention in a functional block diagram, al- 
though all functions, except possibly the first, may be 
carried out by a programmed digital computer. The 
signal S(t a ,t) is first compressed by matched filters 12. ^ 
This achieves range and azimuth correlation in a con- 
ventional manner. The combined operation of these 
matched filters is represented by the first mathematical 
operation 12' in FIG. 4 of correlating the input signal 
S(t a ,t) with a matched filter function f)<t a ,t) to process 3Q 
the input signal at chirp center frequency ftf- 

The outputs e/taojta) from the matched filters are 
then processed through an area filter 14 represented in 
FIG. 4 by the aperture function 14' where {t ao -~ took)- 
A ok, (t 0 —t 0 k)Aok represents a small rectangular area in 35 
the image to be processed that is a square of length \ 0 k 
centered at a point defined by radar times iaokAok in the 
two-dimensional signal S(t a0 ,t 0 ). The result ey*(t a0 , t 0 ) 
for one small area at one wavenumber, ki, is then decor- 
related in block 16 with the inverse matched filter func- 40 
tion for the input signal, and the process is repeated for 
all of the other wavenumbers k2 . . . k n corresponding to 
the frequencies of the matched filters. That decorrela- 
tion for each wave number is represented in FIG. 4 by 
the third mathematical operation 16' where the inverse 45 
matched filter function is represented by f _i . 

Since the procedures of the second and third mathe- 
matical operations (indicated in FIG. 4) are repeated (in 
the apparatus represented by the blocks 14 and 16 in 
FIG. 3) for the same small area at different frequencies 50 
fj, f2 . . . f n , it is necessary for the results to be stored. A 
memory 18 is indicated in FIG. 3 for this purpose. 
There is no corresponding memory indicated in FIG. 4, 
it being understood. 

The next operation represented in FIG. 3 by a symbol 55 
20 is to multiply a pair of signals from the preceding 
operation, such as the signal S/yt(t a <?,t 0 ) for the line k\ 
(frequency fi) with the corresponding signal for the line 
k2, and then multiply the signal for the line ki with the 
signal for the line k3, etc., until every signal has been 60 
multiplied by the signal from line ki to form the product 
of signal pairs Sjk, §/+«,*. This operation of multiplying 
signals in pairs is represented in FIG. 4 by a similar 
multiplication symbol 20'. The products are then trans- 
formed by a Fast Fourier Transform (FFT) 22 and 65 
processed through a squarer 24 shown in FIG. 3. These 
two mathematical operations are similarly represented 
in FIG. 4 by functional blocks 22 and 24. The output of 


the squarer is the spectrum of the input signal corre- 
sponding to one area. 

The corresponding results for each pair of frequen- 
cies (i.e., each ak) of the different images (areas) are 
summed in block 26 of FIG. 3 to increase the signal to 
noise ratio of the system. That is represented by the 
same summing function block 26 in the flow chart of 
FIG. 4. In other words, the block 26 labeled — > sums the 
squares of the products Sy*, Sy +/ |,& for different areas. 

The last operation 28' labeled fo SELECTOR in the 
flow chart of FIG. 4 is to determine the magnitude of 
Doppler shift fcof the output signal from the summing 
operation 26 for any one pair of frequencies for different 
images (areas). The shift determined is of the Ak spec- 
tral line of one pair of frequencies from the center of the 
spectrum. From that shift the ocean current velocity is 
calculated. This complete calculation is represented in 
the functional block diagram of FIG. 3 by a block 28 
labeled v c (Ak) COMPUTER. 

From the foregoing, it is evident that the system of 
the present invention utilizes conventional SAR data 
and processes it to provide images for different sections 
of the chirp bandwidth, equivalent to images obtained 
for different frequencies. A common small area for the 
separate image fields is selected, and the image is decon- 
volved to obtain the signal for the different frequencies 
and same area. A product of pairs of signals are Fourier 
transformed and squared. The spectrum thus obtained 
for different images (areas), but for the same pair of 
frequencies, are added for improved signal to noise 
ratio. 

The spectrum for one area image provides a peak 
shifted from the center by an expected amount due to 
the phase velocity v^= =g/Ak of the Bragg scattering 
wave without current, where g is the acceleration due 
to gravity, and Ak is the chirp center frequency differ- 
ence for the paired signals 9 Sjk and Sy+ n ,jfc. Any differ- 
ence in the amount the peak is shifted from the center of 
the spectrum more than the expected amount is due to 
the radial component of current velocity v c (Ak) as ef- 
fecting the ocean wavelength of lir/kk. 

Computing v c (Ak) is carried out, once fo is deter- 
mined, in the following manner. First it is recalled that 
the phase velocity v^of the ocean surface is known to be 
equal to Vg/Ak, and that the total velocity indicated by 
the Doppler shift f b is given by the following equation 

where f r is the radar frequency, C is the speed of light. 
It is also known that V(Ak)— v^+v^Ak). If V(Ak) is 
thus determined from the Doppler shift and v p is calcu- 
lated, the radial component of velocity v c (Ak) of the 
ocean current resulting from geostrophic forces and 
tides is determined. It is thus seen that the computer 28 
of FIG. 3 performs only simple arithmetic operations 
once fo is determined. The same computations can be 
carried out for each of the other pairs of signals which 
define a Ak, and the computed v t (Ak) of each calcula- 
tion compared or averaged with other calculations for 
greater accuracy. 

The velocity computer 28 operates by effectively 
picking out the Doppler shift fo by a procedure that is 
the equivalent of plotting the Doppler shifted output 
signal from the squarer, and preferably by plotting the 
sum of such signals over different areas to improve 
signal to noise ratio, as just noted above. The computer 
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calculates V(Ak) and v p , and then determines ocean 
current from the equation: 

V{Lk)=v p +v^k) 

This may be converted to ocean current depth profile 
by the relationship: 


v c A/c = 4 k 


p. (z)e -4Afcz d2 


where z is the depth coordinate and p,(z) is the current 
dependence on depth. 

Although particular embodiments of the invention 15 
have been described and illustrated herein, it is recog- 
nized that variations and equivalents may readily occur 
to those skilled in the art. Consequently, it is intended 
that the claims be interpreted to cover such variations 
and equivalents. 20 

What is claimed is: 

1. A method for measurement of ocean current using . 
a chirped synthetic aperture radar which transmits peri- 
odic pulses of predetermined chirp bandwidth toward 
the ocean at an angle Q that is the complement of the 25 
Bragg angle 9, and receives return signals, S(t a ,t), com- 
prising the steps of 

compressing said signals S(t fl ,t) to provide signals 
e/t G 0 ,t 0 ) for different sections of said bandwidth, ^ 
selecting a small common area for separate image 
fields at the different sections of said bandwidth in 
the signals c/tao^o) to produce image signals ejk- 
(t ao,to) for the common area, 
deconvolving the signals ZjkiUoAo) to produce differ- 35 
ent signals for the different frequencies of said 
bandwidth, but of the same common area, 
producing a signal Sjk'Sj+ n ,k that is a product of pairs 
of signals from the previous step, 

Fourier transforming and squaring the signal thus 40 
obtained from the preceding step to produce a 
spectrum, 

determining the Doppler shift fe of the spectrum 
peak from the center of said spectrum, and 
comparing said Doppler shift fc to the Doppler shift 45 
expected due to the phase velocity v p of the Bragg 
scattering wave to obtain ocean current velocity 
v c (Ak) from the difference between the expected 


Doppler phase shift and the determined Doppler 
phase shift. 

2. A method as defined by claim 1 wherein the entire 
process is repeated, each time for a different common 
area, but with the same pairs of frequencies in the signal 
Sjk-Sj+n,k> and the spectra thus obtained for different 
areas using the same pairs of frequencies are added to 
provide improved signal to noise ratio for the spectrum 
used to determine the Doppler shift fc. 

3. Apparatus for measurement of ocean current using 
a chirped synthetic aperture radar which transmits peri- 
odic pulses of predetermined chirp bandwidth toward 
the ocean at an angle ft that is the complement of the 
Bragg angle 0, and receives return signals S(t a ,t), com- 
prising 

means for compressing said signals S(t fl ,t) to provide 
image signals e/Wt 0 ) for different sections of said 
bandwidth, 

means for selecting a common area for separate image 
fields at the different sections of said bandwidth in 
the signals e/t fl 0 ,t 0 ) to produce image signals ejk- 
(t ao ,to) for the common area, 

means for deconvolving the signals e/*(t fl 0 ,t 0 ) to pro- 
duce different signals for the different frequencies 
of said bandwidth, but of the same common area, 

means for producing a signal Sjk • w that is a prod- 
uct of said pair of signals from said deconvolving 
means, 

means for Fourier transforming and squaring the 
signal thus obtained from said means for producing 
a spectrum from said signal %jic'Sj+ n ,k, 

means for determining the Doppler shift fo of the 
spectrum peak from the center of said spectrum, 
and 

means for comparing said Doppler shift fo to the 
expected Doppler shift due to the phase velocity 
v p of the Bragg scattering wave to obtain ocean 
current velocity v c (Ak) from the difference be- 
tween the expected Doppler phase shift and the 
determined Doppler phase shift. 

4. Apparatus as defined in claim 3 wherein said appa- 
ratus repeats operation, each time for a different com- 
mon area, but with the same pairs of frequencies in the 
signal Sjk'Sj+ n ,kf and includes means for adding the spec- 
tra thus obtained for different areas using the same pair 
of frequencies to provide improved signal to noise ratio 
for the spectrum used to select the Doppler shift f d- 
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